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Abstract 

We study coherent enhancement of Coulomb excitation of high energy particles in crystals. 
We develop multiple scattering theory description of coherent excitation which consistently 
incorporates both the specific resonant properties of particle-crystal interactions and the 
final/initial state interaction effects typical of the diffractive scattering. Possible applications 
to observation of induced radiative neutrino transitions are discussed. 



^PHieeQAIX.SP.KFA- JUELICH.DE 
2 ZOLLER@HERON.ITEP.RU 



1 



1 



In a two- level quantum system (the ground state |0) and the excited state |1)) under peri- 
odical perturbation V sin vt with the frequency v equal to the level splitting v\§ = E\ — Eq 
there develop quantum beats with the oscillation frequency uj = (1| V|0). If the perturbation 
V is weak, then for ujt <C 1 the |0) — ► |1) the transition probability P\o(t) increases rapidly 
with the time 

P W (t) OCLUH 2 . (1) 

In particle physics, examples of such rare processes are the weak radiative transitions 
of hyperons, N'j — > Y, beyond-the-standard-model-decays like \i — > ej and radiative (mag- 
netic) neutrino conversion Uij — > v%. For instance, if one could subject hyperons to a high 
frequency field, v ~ my — tjin, the rates of rare decays can be enhanced substantially. The 
high monochromaticity is an evident condition to sustain the growth ([]]) over large time 
scale. 

Okorokov 0,0, P| was the first to suggest that all of the above requirements are met 
best in Coulomb interaction of a high-energy particle propagating in a crystal along the 
crystallographic axis. Here the role of "time" is played by the crystal thickness L. In [|TJ 
the resonant excitation of atoms and nuclei by a periodical Coulomb field of a crystal was 
predicted. In the first observation of the resonant transition of He + (n = 1) to He + (n = 4) 
in Ag crystal was reported. For a propagation through a crystal with the lattice spacing d, 
the frequency v = 2nv/d, where v is the velocity of the atom. For ultrarelativistic particles 
v is enhanced due to the Doppler shift, v = 2iT-fv/d, where 7 is the Lorentz factor. It is 
the Lorentz factor which can boost v to the hundreds MeV's range. Since [|J the Okorokov 
effect has been studied extensively both experimentally ||] and theoretically ||. 

In early works on the subject the Coulomb field of a crystal was evaluated in the 
Weizsaker- Williams approximation [El, [^] and then applied to the calculation of the tran- 
sition amplitude in the plane wave Born approximation. The — > Y transitions in the 
Coulomb field of the nucleus were discussed by Pomeranchuck and Shmushkevich|]] . The 
P7 — > S + excitation in crystals based on the approach [|], and was considered recently 
in and it was claimed that the law Pt, p {N) oc iV 2 , holds up to the crystal thicknesses 
N = L/d~ 10 7 §. 

However, the plane wave Born approximation is not self-contained, it does not allow to 
assess what is the upper limit on iV and, as a matter of fact, it grossly overestimates the 
enhancement factor. Our point is that the coherency of jry — > S + transitions depends on 
the initial state interactions (ISI) of the proton and final state interactions (FSI) of the 
hyperon. It also depends on thermal vibrations of the atoms in a crystal. 

The purpose of this communication is to derive the Okorokov effect for the Pomeranchuck- 
Shmushkevich processes directly from the multiple scattering (MS) theory. We find a dra- 
matic impact of the ISI and FSI effects on p — > S + -transition resulting in the very slow, 
~ log N, rise of the amplitude rather than ~ N. This result is quite general and is equally 
applicable to the excitation of ultrarelativistic nuclei and ions passing through the crystals. 

On the contrary, the magnetic conversion of the neutrino in crystals as shown to have 
a pronounced resonant structure and we comment on possible implications for the future 
laboratory investigations of the neutrino electrodynamics. 



2 



2 



Consider Coulomb interaction of high-energy particle a moving along the < 001 > axis in 
monatomic crystal. The interatomic distances are large, d ~ 5 A, compared to the Thomas- 
Fermi screening radius cltfi a TF — H~ l — (rrieaZ 1 ^ 3 )^ 1 ~ 0.5Z -1 / 3 A, where Z is the atomic 
number and a = 1/137. The relevant impact parameters, b e ff, satisfy b e ff ~ <C d 
and the widely used one-chain approximation is applicable. At the same time, b e ff is much 
larger than the nuclear radius Ra ~ 10~ 5 Z 1 ' 3 A and the effects of strong interactions can be 
safely neglected. 

The amplitude jF ba of coherent transition a — > b on a chain of N identical atoms is 

T ha = (V(r 1 ,...,r N )\S ba \*(ri,...,r N )). (2) 

Here \l/(ri, rjv) is the full ground-state wave function of the iV-atomic chain and r-j = 
(hi, Zi) are positions of atoms. The S-matrix, S ba , admits the expansion into multiple scat- 
tering series, jF 6a = ^2 n F^\ and describes the a-to-b transition as well as both the initial 
and final state interactions of both particles a and b with crystal. The high momentum 
projectile propagates along straight-line trajectories at fixed impact parameter b. Then 
one can use Gribov's dispersion integral representation for the coupled channel n-fold scat- 
tering amplitude F^\h) in the basis of physical states \a), \b) |TD|. This allows one 
to consider the a — > b transition as a sequence of n both diagonal and off-diagonal transi- 
tions | a) — > \i) — > ...\k) — ■> |j) — > \b) ordered along the beam direction. Each off-diagonal 



transition, involves a longitudinal momentum transfer 1 10 



where Am^ = m| — mf , E is the projectile energy and mj is the particle j mass. Associated 
with the Kji is the phase factor exp[in b jZ u + iK,jkZ u -i + ... +iKi a zi], in the amplitude F b ^\b): 

N\ ^ r L 



F£\b) = H)- 1 — — I dz n [d 2 b n f bj (b-b n )eMi Kbj z n ] (4) 

v iV n r i,...,k,j J ° J 

J dz n _ x J Gpb n _i/ i/c (b - b n _i) exp[zK jfc z n _i]... 
... J dzi y<i 2 bi/ ia (b-bi)exp[2fi; ia zi]x(ri,...,r n ), 

where 

X(ri,...,r n ) = J d 3 r n+1 ...d 3 r N \^(r r N )\ 2 . (5) 
For our purposes it is sufficient to use the uncorrected wave function 



jV 



|*(r 1 ,...,r JV )| 2 = n^("<)» ( 6 ) 



i=l 



where the normal coordinates are defined by = (z — l)d + Uj, % = 1, N, d = (0, 0, d) 
and u, = (bi,u zi ). 

For the rare processes of interest 

fij < fujjj ,i ^ J ( 7 ) 
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and we keep only the lowest order terms in powers of the off-diagonal transitions. Then, 



N N m-1 N-m 

^a(b) = £ F£>(b) = E E E h aa (brh bb (brh ba (b)St ] (K ba ) , (8) 

n=l m=l ni=0 n2=0 



where 



h ba (h) = J d 2 hif ba (h - bj)p(bi) 



P(bi 



du~i(b( Ui 



S < £ n) (K ba )= I au zm (p{u m )exp{iK ba u zm , 
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du 7 



2tt(u 2 ) 



cxp 



exp 



3b? 



2(u 2 ) 



(9) 
(10) 



6 



exp [iK ba (m - l)d] . (11) 



For the sake of simplicity, we have assumed the Gaussian form of 0(u«). In ([10|) (u 2 ) is 
the mean squared amplitude of thermal vibrations of the lattice . Note that the relevant 
longitudinal momenta are as follows Kji ~ d~ l <C (u 2 ) -1 / 2 . The summation over m in eq.(H) 
results in 



F b a(b) = h ba (b)exp 
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(12) 



x 



[1 + th bb {b)] N - [1 + ih aa {b)} N exp [iKbgNd] 
[1 + ih bb (b)] - [1 + ih aa (h)} exp [iK ba d] 



In the practically interesting case like jry — > S + and/or e7 — > we have / aa = and the 
amplitube jF fca reads 



^(b) = /l6o(b) [1 + l/loaCb)]"^ 1 S L (K ba ) 

where Sl(k) is the longitudinal form factor 



SL(nba) = exp 

Evidently, for the resonant values of K ba , 

27m 

K-ba 



2 / 2\ 



sin(/t ba iVrf/2) 
sin(/t 6a d/2) 



n = 0,±1,±2. 



we have Sl oc iV. The resonance energy E n , 



Am 2 ba d 



Ann 



can be lowered going to n 1. However, there is an upper limit on n 



1 (3d 



2 \ 1/2 



2tt V(u 2 ), 



5-10 



(13) 



(14) 



(15) 



(16) 



;i7) 



which comes from the Debye- Waller factor in Sl(k), eq. (|I~4|) . 

Above we focused on the S-matrix at a fixed impact parameter b. The experimentally 
observable quantity is the differential cross section da ba /dqj_ = \T ba (q±)\ 2 , where 



r d 2 h 

%a(q±) = S L (K ba ) / — - h ba (b) [1 + z/i^b)]^ -1 exp[iq ± b] 
J Zn 



18) 



and q_i_ - is the two-dimensional vector of the transverse momentum (|qj_| = q±). 

As we shall see below, the dominant contribution to (pT8|) comes from b e ff ~ fi~ 1 N e ff, 
where 

N eff = log(ZaN) > 1 . 

Then, because (u 2 ) ~ /i~ 2 

h aa (h) = J d 2 bj aa (b ~ b«)p(b l ) ~ f aa (b) . (19) 

The amplitude f aa of elastic scattering of a charged particle a in a screened Coulomb field 
of an isolated atom is purely real if q\ <C /J,p, where p is the projectile momentum: 

f aa (b) = ZaK (fib) . (20) 

Here we neglect the contributions of anomalous magnetic moments to the small-angle elastic 
scattering. 

Let the resonance condition ( |T5| ) be satisfyed so that Sl(k) oc N. Our point is that 
because of the distortion factor [1 + ih aa (b)] N in the integrand of eq . (|l"8D the law T ba oc N 
does not hold. The evaluation of the distortion effects for the p7 — > S + transition proceeds 
as follows. 

For the phenomenological transition matrix element P|, [Tl"| 

M(pj -> S + ) = iuy, (/iE P + 75^s P ) cr^e^Up (21) 

one readily finds 

/sp(b) = [/ispcr [b x n] + id^ p crh] b' 1 Z^/apK^b) . (22) 

Here n - is a unit vector along the projectile direction, cr is the Pauli spin vector and Ki(x) 
is a modified Bessel function. 

Then, the evaluation of the full transition amplitude reads 

^E P (q±) = [ifJ'SpCT [q± x n] + d^ p aq ± } q]} /iS £(re ba ) (23) 
xZ^ J bdb JxiqxfyKxdib)) [1 + iZaK (^b)} N ' 1 

and the steepest descent from the saddle-point at fib = N e ff — in/2 yields for q± < f^N~^ 
T Sp (g ± ) oc i (N eff - m/2) J, (q x ^ 1 {N ef} - m/2)) . (24) 



where J\{x) is the Bessel function. Consequently, as soon as NaZ ^> 1, which holds for all 
practical purposes, even under the resonance condition (|15|) one only has the logarithmic 
growth 

T Sp oc \og(ZaN) . (25) 

Although the relevant impact parameters rise, 6 e // ~ [i~ l N e ff, for all the practical pur- 
poses b e ff <C d so that the one-chain approximation holds. The total cross-section, a = 
J dqj_\TY, p (qj_)\ 2 appeares to be independent of N. 

We conclude that ISI and FSI effects completely destroy the resonant enhancement effect. 
Similar phenomena are well known in the high-energy diffractive scattering of hadrons, 
though for purely imaginary amplitudes. For the screening of real elastic amplitudes in the 
Coulomb potential scattering see e.g. Ul2"|. 
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Several proposals [13) of experiments on the neutrino conversion induced by the neutrino 
magnetic moment in external electromagnetic field 



1S 2 



(26) 



are under discussion [14]. Hereafter we are dealing with the non-diagonal magnetic (/i2i) and 
electric (d 2 i) dipole moments which exist in general for both Dirac and Majorana neutrinos 



PI 



The experimental bounds on /i 2 i are such fl4| that the FSI/ISI effects are negligible 
for all practical purposes. Then, in terms of the matrix element (pT|) the differential cross 
section of the process (p6|) reads 



da- 



21 



dq 2 L 



aZ z b: 



2l2 
21 



(ql + /i 2 ) 2 



(27) 



where b\ x 



|/^2i | 2 + |^2i| 2 for Dirac neutrinos and b\\ = 4 [(Im/i 2 i) 2 + (Rec^i) 2 ] for Majorana 



neutrinos. The transverse form factor Sx(qj_) equals 



exp 



2 / 2\ 



(21 



The observation of the Uij —>■ z/ 2 -transition which would look like a resonance at some 
neutrino energy E ~ E n (see eq. fllB"!) ) would enable one to evaluate the mass difference of 
the neutrinos of different species and estimate the neutrino transition magnetic moment. 

For these purposes it is useful to represent the resonant part of the differential cross 
section ( PT| ) as a sum of Lorentz curves. Then, the neutrino conversion rate, R c = 0"2i/gP, 
equals 



R r ~ N 2 



blaZ* 



d? 



1 + |/i 2 (u 2 ) 



E 



(r B /2) s 



The width of n-th resonance is 



(r n /2) 2 + (E-E n 



V3 Ami d 



exp 



47r 2 n 2 (u 2 ) 
~~ 3 dz~ 



(29) 



(30) 



2 n 2 n 2 N 

Then, the " disappearance type" experiment [0 with the CERN 30 GeV z/-beam sup- 
plemented with the crystal target would allow to investigate the region of Am 2 x up to 
~ (lMeV) 2 . If we accept that the departure of the neutrino from the crystallographic 
axis which does destroy the coherence is due to the transverse momenta in the range from 
q± ~ d' 1 up to q± ~ fi, the crystal-target experiment could hope to detect transition 
neutrino magnetic moments larger than 



A*21 



:io- 



10~ 8 )^ 



B, 



(31) 



where fiB is Bohr magneton. The efficiency is assumed to be ~ 10 -4 |T3j. It should 
be emphasized that in the v — crystal coherent interaction typical q 2 vary in the range 
~ (lOKeV) 2 . For comparison, presently available bounds on fi 2 i were obtained for virtual 
photons with typical q 2 ~ (100 MeV) 2 . A nontrivial q 2 dependence, resulting in a strong 
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low-momentum enhancement of a neutrino magnetic form factor, discussed in [T3) makes 
such an experiment topical. 
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